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A noise source model, consisting of a pulse sequence at random times with memory, is 
presented. By varying the memory we can obtain variable randomness of the stochastic 
process. The delay time between pulses, i. e. the noise memory, produces different kinds 
of correlated noise ranging from white noise, without delay, to quasi-periodical process, 
with delay close to the average period of the pulses. The spectral density is calculated. 
This type of noise could be useful to describe physical and biological systems where 
some delay is present. In particular it could be useful in population dynamics. A simple 
dynamical model for epidemiological infection with this noise source is presented. We 
find that the time behavior of the illness depends on the noise parameters. Specifically 
the amplitude and the memory of the noise affect the number of infected people. 

Keywords: Statistical mechanics; population dynamics; noise induced effects. 

1. Introduction 

Many studies have reported occurrence of fluctuations and noise in biological sys- 
tems and the role played by the noise, always present in natural systems [1-3]. 
The stochastic nature of living systems and the related noise-induced eS'ects ranges 
from bio-informatics [4,5], population dynamics [6-9], to virus dynamics and epi- 
demics [10, 11]. Noise existing in biological systems is due to environmental fluctu- 
ations and it is usually taken into account as a multiplicative white noise [12]. The 
white noise is usually deflned as a continuous stochastic process with zero mean and 
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correlation function (^^t) — C6(t), where C is the noise intensity. However when 
we reproduce this noise in discrete time with step r for computer simulations, we 
consider rectangular pulses with width r and some distribution of probabilities for 
their heights. In this case the randomness is in the amplitude of noise pulses. We 
call this noise a discrete white noise source. In this paper we will consider both 
continuous and discrete noise sources. Starting from Markovian white noise we 
will introduce some memory obtaining a quasi-periodical process, which is strongly 
NonMarkovian . 

In some population dynamics models, the noise must be positive at all times and 
often the randomness is due to stochastic natural events, that can be represented as 
pulses occurring at random times. We consider therefore a white noise source as a 
series of positive pulses at random times. To obtain zero average, if it is necessary, 
we can subtract the positive mean value. 

In our model the discrete white noise at i-th step can be expressed as 

T 

where //r is the amplitude of pulses, which we set constant for simplicity, rji is 
a random value, distributed homogeneously in the interval [0, 1], 'd{x) is the theta 
function > 0) = 1, t9(x < 0) = 0), v is the probability per unit time to have a 
pulse and vt is the probability per one step to have a pulse. To obtain a time series 
of noise we extract random values rji at every time step r. The average distance 
between pulses C, can be considered as an effective period (C) = T = Xjv (with 
T ^ t). The noise average is (f) = f /T, and the correlation function is different 
from zero only during the step r. From formal point of view this noise is therefore 
white. 

The continuous white noise is obtained in the limit r — > 0. This noise is usually 
expressed as a sequence of i5-shape pulses at random times tj [13] 

3 

where 

=io + Ci + -- + 0, (3) 

with ^0 the initial time, and Q the random time distances between neighboring 
pulses. The probability distribution is w(C) = i^exp(— i/C), where C is distributed in 
the interval [0, oo]. Of course the effective period T ~ {Q and the mean value (5) are 
the same as for discrete noise, and here we consider them as constant parameters. 
We will use both continuous 10) and discrete |^ descriptions of the noise. 

2. Noise vifith memory 

For some biological applications should be interesting to consider a noise source 
with some delay C.o after each pulse. The system has memory during time Co and 
the next pulse is forbidden during this time. This process is suitable to obtain 
noise sources with varying degree of randomness, ranging from white noise, as we 
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defined it above ((^o = 0), to quasi-periodical process (Co — T). So the probability 
distribution of random time distances Cj is 

it;(() =pef('^°~«), (4) 

where C is distributed in the interval [Co , oo] , p is the probability per unit time to 
have a pulse. If Co =0, then p = v. The first moment of C is therefore easily derived 
as 

(C)-Co + -. (5) 
P 

For fixed effective period (C), the probability density p increases as memory Co 
increases. This probability is equal to 1/r, when we have periodical process in 
discrete time (Co — > (C) ~ ''")• For continuous time, the probability of pulse per unit 
time p — > CX3, when Co ^ T. 

For the discrete process the standard deviation of the noise is the weighted sum 
of this quantity during delay (which is zero) and that after delay a"^ , with the related 
probabilities. All parameters after delay we will mark by tilde. So we obtain 

a^ = |xO+^xa^. (6) 
For the variance after delay we have 

^^ = (l^)-(l>^=({)^T-(4r)', (7) 

where (^) and (f^) are taken for the time after delay. Substituting Q into © and 
taking into account (jSJ we obtain 

So we have an increasing as the memory of the system decreases. The standard 
deviation achieves its maximum value ~ ((T — t) / [T"^ t)) f'^ ~ P /tT for Co — 0. 
When the memory is maximum C^f) ^ T — t, the process is not random, so cr^ = 0. 

For the continuous noise process it is useful to calculate the spectral density. 
The characteristic function for C is the Fourier transform of the probability density 

OO 

g-K+«-Crf^ = i- . (9) 

p-iuj 

Co 

By using the relation between the spectral density of the process ^ and the charac- 
teristic function of C, derived in ref. [13] 

2(1-|9(^)P) 

we get 
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2^2(1 - Co)' 



(11) 



2-2 cos(tjCo) + 2(1 - Co)w sin(t^Co) + ^^(i _ ^p)2 ' 
where we use T =< ^ 1- The power spectrum of the noise ^(t) is shown in 
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Fig 1. Power spectrum of the process ^(t) for three values of the memory (^q: curves 1, 2 and 3, 
correspond respectively to = 0.99,0.3,0.01. 

Fig. n for three values of memory parameter, namely Co = 0.99,0.3,0.01. We see 
that at Co = 0.99 the spectrum is a series of J-function-like pulses with harmonics at 
uJn — 27m. For Co — 0.01 the spectrum density is almost constant, and corresponds 
to white noise. For all cases, the spectrum density is quasi-periodical with distance 
between maxima equal to: ujn — LOn-i — 27r/Co. 



3. Some possible applications 

The pulse noise model can be useful to describe some different physical and bio- 
logical systems, where random processes with delay are present. It can be useful 
to describe and to understand the behavior of an epidemiological infection inside 
a population with a big number of individuals. In particular the number of con- 
tacts between ill and healthy people can be modelled as our noise with memory, 
where the delay describes the incubation period of the illness. So we can compare 
development of illnesses with different incubation periods. 

Another possible application is the neuronal response to an electrochemical 
pulse. In fact the response time of a neuron, activated by a pulse, can depend 
on temperature and other parameters. Moreover our pulse noise model allows to 
regulate the periodicity, and this should play an important role in the neural net- 
work field. 

The pulse noise source can also represent an alternative way to describe popula- 
tion dynamics of interacting species, previously modelled using stochastic resonance 
phenomenon [8] . This should be very useful when a periodicity, connected with en- 
vironmental parameters, is present. 
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4. Simple dynamics in the presence of noise with memory 

We start to analyze the role of the noise with memory considering a simple system 
described by the equation 

+ (12) 

where ^(i) is the the noise with memory, a is a negative constant, and I{t) can 
be a species density in population dynamics or the number of infected people in 
epidemics. We interpret this stochastic differential equation in the Ito sense. This 
equation has been previously investigated in refs. [14, 15]. 
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Fig 2. Pulse noise ^(t) and corresponding time series of the number of infected people lit) vs 

time, for two values of the delay C,q. The values of the parameters are: T = {Q = lOr, r = 10~^, 

/ = 3 . 10-2, 7(0) = 0.5, and a = -2. (a) white noise source, with = 0; (b) correlated noise 
source, with = St. 



The formal solution of Eq. ((T^ is 

/(t) = l{0)exp 



at 



(13) 



By numerical integration of Eq. p2(l we find that the noise amplitude parameter 
/ has a critical value fc- For f > fc &n instability occurs and the illness exhibits 
a divergent behavior. By setting constant the value of the average period T = 
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(C) = lOr, time series of the infection density are obtained for different values of 
the delay which represents the time before a new pulse can occur. We report 
the results in Fig. |21 together with the pulse noise source. The parameter setting 
is: /(O) = 0.5, / = 3 • 10"^^ a = -2 and r = 10"^. In Fig. El we show the same 
quantities £,{t) and I{t), for the same parameter setting with / = 10~^. In both 
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Fig 3. Pulse noise ^(t) and corresponding time series of the number of infected people I{t) vs 
time, for two values of the delay i^o- Here / = 10~^, and all other parameter values are the same 
of Fig|^ (a) white noise source, with (^o = 0; (b) correlated noise source, with (^o = St. 





figures an exponential behavior of I(t) appears with fluctuations more pronounced, 
when the noise source is white (^o = 0) (see Figs.EK and (St). For correlated noise 
we have a quasi-periodical process. The exponential behavior is decreasing (see 
Fig. [SJa) or increasing (see Fig. ^) depending on the value of parameter / with 
respect to the critical value fc- We observe clearly that the noise memory affects 
directly the fluctuations of the number of infected people I{t). Using noise with 
zero memory we find that the critical value is /c = 2 • 10~^. The noise memory also 
affects the exact value of fc- In particular, when / is near the critical value with zero 
memory, the time behavior of the average number of infected people changes from 
increasing to decreasing exponential behavior, depending on the value of the noise 
memory. In Fig. 01 we report the time behavior of the average ensemble < I{t) > for 
different values of noise parameters / and Cq. Namely for (a) Co = Oi / = 2.01 • 10^^, 
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(b) Co = 0, / = /c = 2.00 • 10-2, and (c) Co = 8t, / = 2.01 ■ lO^^. The average 
is performed on lO** numerical experiments. We see that for / — 2.01 • 10"^ ~ /c 
(critical value for white noise with zero memory), by increasing the noise memory 
from zero to Co = 8r, an instability-stability transition occurs. The average number 
of infected people changes from increasing ((a) in Fig.^ to decreasing ((b) in Fig.^ 
exponential behavior. 
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Fig 4. Time behavior of the ensemble average of the number of infected people I{t) series, for 
different values both of the delay and the noise amplitude /. All the other parameter values 
are the same of Fig|21 (a) white noise source (Co = 0) with / = 2.01 ■ 10-2; (b) white noise source 
(Co = 0) with / = 2.00 • 10-2; (c) correlated noise source with Co = 8t and / = 2.01 • 10-^. The 
number of numerical realizations is 10*. 

5. Conclusions 

The noise with memory described above is a very useful tool for modelling quasi- 
periodical processes in nature. The parameter of periodicity allows to distinguish 
deterministic and random processes, and to adjust model to a real process accord- 
ingly. When this noise is included into the equation as a multiplicative noise, we 
can see qualitatively different behavior of the system in dependence on the noise 
amplitude and the periodicity parameter of the noise source. The role of the noise 
memory is to induce an enhancement of fluctuations of the number of infected peo- 
ple I{t) and to shift the critical value. When the noise parameter / is near the 
critical value with zero memory, a stability-instability transition occurs. For epi- 
demic dynamics we find that, for pulse amplitude parameter less than the critical 
value (/ < /c), we obtain quick recovering. Moreover we find that the role of the 
noise memory, that is the incubation period, is to reduce the amplitude of the fluc- 
tuations of the infected people. This causes a more regular time behavior of the 
illness. 

Finally concerning population dynamics of some species, this noise can represent 
situations in which the species density growth is a quick process, and the decreasing 
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process is continuous in time. The delay Co should represent the time taken by the 

individuals of the species to become adult. So wc could estimate the dependence of 
population dynamics on this parameter, with all other ones fixed. 
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